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Abstract. Proteus mirabilis are bacteria that make strikingly regular spatial-temporal patterns 
on agar surfaces. In this paper we investigate a mathematical model that has been shown to 
display these structures when solved numerically. The model consists of an ordinary differential 
equation coupled with a partial differential equation involving a first-order hyperbolic aging 
term together with nonlinear degenerate diffusion. The system is shown to admit global weak 
solutions. 
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1. Introduction 

Bacteria of the species Proteus mirabilis are ubiquitous throughout nature. In human beings, 
Proteus mirabilis is found as part of the normal flora of the gut. Its main pathological role is 
in infections of the urinary tract, but it can also cause wound infections and septicaemia. Even 
though most of the human urinary tract infections are due to the bacterium Escherichia coli, 
urinary tract infections due to Proteus mirabilis are also well-documented. It commonly invades 
the urinary tract when the normal function of the tract is disturbed by instrumentation such 
as catheterization. Once attached to urinary tract, Proteus mirabilis infects the kidney more 
commonly than Escherichia coli and characteristically leads to urinary stones. 

Proteus mirabilis can exist in two distinct morphological and physiological forms known as 
"swimmer" cells and "swarmer" cells, respectively. Broth cultures of Proteus mirabilis consist 
virtually exclusively of mononuclear cells (swimmers) approximately 1 fim wide with short flagclla. 
Swimmer cells go through a prototypical bacterial cell growth and division cycle. However, when 
inoculated onto agar surfaces, some cells cease septation but continue to grow and produce many 
lateral flagclla to form elongated multi-nucleoid hyperflagellated swarmer cells up to 100 urn in 
length which aggregate in parallel arrays to form motile multicellular "rafts" . The process in 
which dividing cells become swarmers is called "differentiation" and occurs only above a critical 
dividing-cell density. Rafts of swarmer cells are capable of translocation while swimmer cells are 
immobile. The movement of Proteus mirabilis through raft building requires two things, namely 
sufficient maturity in swarmer cells to contribute to raft building and sufficient biomass of mature 
cells to form the rafts. After some time migrating, when the multinuclear swarmer cells approach a 
maximal size, they cease movement and rapidly "dedifferentiate" again into single nucleus swimmer 
cells. This coordinated burst of swarming activity interspersed with a consolidation to the swimmer 
state results in characteristic concentric rings of growth. It is due to these strikingly regular spatial 
and temporal patterns that Proteus mirabilis has attracted attention in the mathematical biology 
literature [SEE [12]. 

The ability to form swarmer cells seems to allow rapid colonization of solid surfaces and the 
establishment of extensive Proteus mirabilis biofilms. Apparently, Proteus mirabilis can also swarm 
over the surfaces of all the major catheter. Swarming may thus play roles in both the initiation of 
catheter associated infections and the subsequent spread of the biofilm over the catheter surface. 
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The processes involved in the evolution of Proteus mirabilis and the formation of regular patterns 
are rather complex. A key ingredient of the mathematical representation is the age dependence of 
swarmer cell behavior. An age and spatially structured model for Proteus mirabilis swarm colony 
development was presented in [5] , and - in slightly modified form - in |12j and [51 [7] . Denoting by 
v = v(t,x) the swimmer cell density in dependence of time t > and spatial position x S (with 
a spatial region f2 C W 1 ) and by u = u(t, a, x) the swarmer cell density which additionally depends 
on age a > 0, the models in [51 [7J Q~2] can be re-cast in the form 

d t u + d a u = div x (D(A(t,x)) 

9 t V = - (l - £(v)) v + / e a/r Li{a)u{t,a,x)da , (t, x) G (0, oo) x Q , (2) 
r Jo 

u(t,0,x) = -£(v(t,x))v(t,x) , {t, x) e (0, oo) x fi , (3) 

w(0, a, x) = u°(a, a;) , v(0, a:) = v°(x) , (a, x) e (0, oo) x O , (4) 
d u u = , (t,a,x) <E (0, oo) x (0, oo) x <9S1 , (5) 

where 

/•OO 

A(t,x) := / e a/r a, x) da , (t, x) e (0, oo) x O . (6) 

The major differences of the models [BJ O and their philosophies are in different choices of the 
functions D, fx, and £. The meaning of the various terms are as follows: A = A(i, x) represents the 
total motile swarmer cell biomass, where ao > is the minimal age of swarmer cells required to 
participate actively in group migration. The exponential comes in since biomass increase during 
swarm development is assumed to occur at the same rate as during the swimmer cell cycle. The 
parameter r is the time it takes a cell to subdivide. 

Equation ([1]) expresses the change in time of swarmer cells of a given age a. Movement of 
Proteus mirabilis occurs if sufficiently many swarmers above the critical age ao group together to 
build a mass above a certain threshold A m j„ > 0. Thus, the diffusivity D depends on A and is 
small (or zero) for A small. For instance, D may be of the form 

D(A) = D max{A - A mm , 0}™- 1 (7) 

as in [5] with D E (0, oo) and m = 2. In [5] also a dependence of D on v and on a memory 
term is included, something we will refrain of taking into account. Note that the exponential 
weighting in A means in ([7J that older cells contribute more to swarming than younger cells. The 
age dependent function [i in (fTJ) is the dedifferentiation modulus, which is higher for older swarmers 
than for younger ones. A typical shape for fi is a narrow hump located around a maximal age 
a-max and zero elsewhere. The limit choice /i(a) = fio <5 a =a max has also been considered in [9]. 

The change in time of the swimmer population is given by equation ([2|) . The population grows 
exponentially with rate 1/r. Some of the swimmer cells cease septation and differentiate with 
rate £,(v)/t into swarmers of age 0. This increase in swarmer cells is reflected by equation ([3]). As 
pointed out in [6] the function £ should be zero for v small. Indeed, the incorporation of a lag phase 
in swarmer cell production triggers the development of a consolidation phase after a swarm phase 
and thus prevents a self-sustaining soliton caused by swarmers that dedifferentiate into swimmers 
immediately differentiating into new swarmers. This lag in the onset of differentiation was observed 
in fTJ] and included in the models in [7J [T^] . The integral term in ^ represents dedifferentiation 
of swarmer cells into swimmer cells. 

The basis of equations CE])-© was presented in [9] and extensions and modifications of these 
equations were proposed in [6l [Jj [12] . In [9] and [H [7] the main focus - besides the modeling 
aspect - were computational results displaying the spatial and temporal patterns of concentric 
rings with equal width. In [7j numerical results were presented examining the necessity of a 
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sharp age of dedifferentiation from swarmer to swimmer cells. All papers [5J [7J [5] use explicit age 
dependence in the evolution of the swarmers. As pointed out in [5] explicit age structure provides a 
mechanism for controlling - at least numerically - the ratio of time spent swarming to time spent in 
consolidation without changing the total cycle time. In [12] a reaction-diffusion model for Proteus 
mirabilis swarm-colony development based on averaging over the age variable was used and results 
on the long time distribution A/v were derived. A model for the periodic swarming of Proteus 
ignoring the age structure from the outset was introduced in [5]. 

For further reading concerning morphology and pathogenicity of Proteus mirabilis and for nu- 
merical results for models of this bacteria we refer to [6j [Jj [9l [T2l [13l [T4J [TTJ [18] and the references 
therein. 

On the other hand, fewer mathematical results seem to be available for models of Proteus 
mirabilis and the only result regarding the mathematical well-posedness of models for Proteus 
mirabilis we are aware of is |10) . Existence and uniqueness of weak solutions to ([I])-© are shown 
in [10] for the case of non-degenerate diffusion including memory. 

The purpose of this paper is to prove an existence result for diffusion coefficients D(A) that may 
degenerate for A = and thus to get closer to the biological reality. Still we cannot handle the 
case where D is given by (J7J) but expect that the outcome of the model with diffusivity 

£>(A) = Do max{A - A min> 0}" 1 - 1 + e' 1 ^ A > 

(to which our result applies if m > 3) for small s > resembles that for e — from a numerical 
viewpoint. However, proving the formation of regular spatio-temporal patterns is beyond the scope 
of this paper. 

The outline of the paper is as follows: In the next section we first establish an existence and 
uniqueness result for the non-degenerate case; that is, when D is bounded below by a positive 
constant. Our method for proving this result is completely different from that in [10] . Section [3] 
then shows how to handle certain degenerate diffusivity. 

Throughout the paper we assume that the minimal age ao required for swarmer cells to partici- 
pate actively in the collective motion is positive. The case ao = turns out to be easier and could 
also be handled with minor modifications. 

2. The Non-Degenerate Case 

Throughout this section we suppose that the diffusivity D satisfies 

DeC 2 ~(R) and D(z) > al Q > for zel, (8) 

where C k ~ (resp. C*~) for k G N\ {0} denotes the set of C fc_1 -smooth functions with a Lipschitz 
continuous (resp. uniformly Lipschitz continuous on bounded subsets) (A; — l)-th derivative. As 
for the differentiation rate £ we assume that 

£eC* 3 -(K) and < £(z) < 1 for zel, (9) 

while the dedifferentiation modulus 

fj, G BC(R) := C(E) n Loo(M) is non-negative . (10) 

Let fl be a bounded and smooth domain in E n . We fix p > n and denote by either the 

space Wp a := Wp a (tt) if 2<r < 1 + 1/p or the subspace of W^ a {0) consisting of those elements 
satisfying homogeneous Neumann boundary conditions if 2a > 1 + For abbreviation we put 
E a := L 1 {R+,W^ B , e a / T da) for a G [Q, 1]. 

In the following we denote by c(T), c(i?), and c(T,R) constants depending increasingly on the 
arguments and that may differ from occurrence to occurrence . 

We first prove an auxiliary result regarding the solvability of ^ . 
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Lemma 2.1. Let T > 0, 2<r G (1 + n/p,2), and assume that v° G W* B and u G C([0, T],E a ) 
are given non-negative functions. Then there exists a unique solution v — v u G C ([0, T],Wp%) 
to ([2|) subject to the initial condition v(0) = V . This solution is non-negative and belongs to 
C([0,T],W* B ) ifueL^TlEr). 

Moreover, if u and u both belong to C([0, T], E a ) and satisfy max {||«(t)||.E . , ||w(i)||.E CT } < R for 
t G [0,T] and some R > 0, then 

\\v u (t) - v u (s)\\ W 2, <c(T,R)\t-s\ , 0<t,s<T (11) 

and 

\\v u {t) - v % {t)\\ wr <c{T,R)\\u-u\\ c{[q ,t],e„), 0<t<T (12) 
for some constant c(T, R) > 0. 

Proof. First note that the regularity of £ and [3J Theorem 4.2] imply 

[v -> 6 C£" , W p 2 £) (13) 

since 2<r > 1+n/p. In addition, (fit))) ensures that the integral term in ([5]) belongs to C([0, L], VF 2 g). 
The existence of a unique non-negative solution v u G C 1 (J, W^g) is now obvious, where either 
J = [0, T] or J = [0, f) with f < T and \\v u {t)\\ W 2„ -> oo as t / f. Next, © and the embedding 
Wp CT <^-> Lqo ensure w„ G L^J, Loo)- Then, taking the gradient with respect to x on both sides of 
© we similarly obtain v u G L OQ (J, W^,). Recalling that pointwise multiplication satisfies 

according to [3j Theorem 4.2] since 2er — 1 > n/p we deduce that 

UMt))v u (t)\\ W 2, < c(\ + \\v u {t)\\ wr ) , teJ . 

^From this we first conclude that v u G L 00 (J,W^ a ) 1 whence J = [0,T], and then £(v)v G 
Loo(L Wl a ) so that (JTTJ) follows by ©. Property (JT^J is implied by dTUJ) and (JT3J>. Finally, if 
u G Li([0,T],Li), then the integral term in © belongs to Li([0, T], W^ B ) due to (TTUD, and we 
readily infer that v u belongs to C([0, T], W* B ). □ 

The solvability of |T|) is based on the following formal observation: Suppose that the function u 
is sufficiently smooth so that the function A = A„, given by ©, leads to a well-defined evolution 
system Ua u (t, s) on L p corresponding to the differential operator 

A u {t)w:= -div x (D(A u (t))V x w) , w G % 2 e . (14) 

Then |T]), (PJ-(l5]) can be re- written as a problem in L p of the form 

d t u + d a u + fx(a)u = -A u (t)u , a > , < t < T , 

u(i,0,-) = -£(«„(<)) « u (t) , 0<i<T, 
w(0,a, •) = u°(a, •) , a > , 

where d u is the corresponding solution to p]). Applying the method of characteristics we derive 
that u is a fixed point of the map $, given by 

( -e-ti^ dr U Au {t,t-a)£(v u (t-a))v u (t-a) , 0<a<t, 
*(«)(«, o) := ^ T (15) 
[ e - /.'-« ^)dr ^ (t,0) u °(a-t) , 0<t<a. 

We now show that this map <& indeed has a fixed point in a suitable space and thus (H}-(0 admits 
a unique solution. More precisely, we have: 
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Theorem 2.2. Suppose ([5|)- (fTu]) and fixp > n, q > 1 and2uj G (1 + n/p, 2). Consider non-negative 
initial values v° G W^g and 

M ° G Si n C\[0, oo), L p ) n L,((0, ao), w^f) n C([0, a ], ^p 2 "" 1 ) 

satisfying the compatibility condition £ > (v°)v° = tu (0, •) in f2. TTien problem (|T|) -((5J) possesses a 
unique non-negative solution (v, u) with 

v G C^flO, oo), W^) n C([0, oo), W 2 v3 ) , u G C([0, oo), £„) n L^ttO, oo), Si) 

/or any 77 G (0, 1) and smc/i that u satisfies 

d t u(t, •) , d a u(t, •) g c([o, t],L p ) n c((<, oo), l p ) , 

d t u(-,a) , <9 a it(-, a) G C([0, a), L p ) n C([a, 00), L p ) 

for all t, a > and solves ([T]) in L p /or t =/= a. 

Proof. Given rj G (0, 1) we fix numbers c, and p such that 

1 + n/p < 2d < 2uo < 2a < 2g < 2 

and rj < a and choose k G (0, minjo" — 1?, 1/q'}), where q' is the dual exponent of q. Note that we 
may assume without loss of generality that quj < 1 by making q smaller if necessary. Let Co be the 
norm of the natural injection W p g c — > W p a B and let R > be such that 

coe 1/T \\u°\\ Ee + ||« |U,((o,«o) )W? -) < R ■ (16) 

For T G (0,1) we denote by Vt the space consisting of all non-negative u G C([0, T], E a ) such 
that ||u(£)|| E(T < R + 1 and ||A„(i) - A„(s)||v^2« < |t - s| K for < t, s < T, where A„ is given 
by ([6]). Then, given any u G Vt, it follows that the operator —A u (t) defined in (fl4|) is for each 
£ G [0, T] the generator of a positive analytic semigroup on L p (e.g. [TJHB]). Moreover, due to the 
embedding W p CT W p § we have 

\\A u {t)~A u {s)\\ c[W 2^ Lp) <c(R)\t~s\ K , 0<t,s<T, 

and 

\\A u (t) - A ri (t)\\c(wl B ,L p ) < c(R)\\u-u\\ Vt , 0<t<T, u,SgV t , (17) 
with the notation 

\\u\\v T := sup ||tt(i)||.E„ for u G Vr ■ 
te[o,T] 

Therefore, invoking Corollary II. 4. 4. 2, Lemma II. 5. 1.3, Lemma II. 5. 1.4, Equation (II. 5. 3. 8), and 
Section II. 6. 4 in [1] and using standard interpolation results on Sobolev spaces with boundary 
conditions we derive that, for any u G Vt, there exists a unique positive evolution system Ua u (t, s), 
< s < t < T on L p such that 

\\UA u (t,s)\\ C (w* % ) + (t-sr- a \\U Au (t,s)\\ c{w2 , < c(R) (18) 

for0<s<t<T, 0<a</3<7<l with 2/3, 2j ^ 1 + 1/p, and 

||f/ Al ,(t,r)-C/^( S ,r)|| } < c(i?)(i- S )^ (19) 

for 0<r<s<<<T, 0</3<7<l with 2/3, 27 7^ 1 + 1/p. In addition, if u is another function 
in Vt, we have 

\\U Au {t,s)-U Aa {t,s)\\ aw2a W 2 0) < c{R)(t-s) a - fJ \\u-u\\ VT (20) 
for < s < t < T, < a, (3 < 1 with a ^ 0, (3 ^ 1, 2a, 2/3 ^ 1 + 1/p. 
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Since — A x subject to homogeneous Neumann conditions on the boundary generates a contrac- 
tion semigroup on W£ g according to [H Corollary V.2.1.4] it follows from (j2"0|) that 

\\ u A.fo*)h(w&,w% il ) ^ c(R) (t ~ sy~° + Co (21) 
for < s < t < T. Also note that 

l|0A.MIU(z r ) < 1 , r6(l,oo], 0<s<t<T. (22) 

Defining $ by (fT5|) we now claim that $ : Vt — > Vt is a contraction provided T = T(R) G (0, 1) 
is chosen sufficiently small. To prove this we fix u G Vt and observe that v u G C 1 ([0, T], is 
well-defined due to Lemma 12. f I Furthermore, Lemma 12.11 and (fT5|) entail 



U(vu(t))v u (t)\\ W 2« < c(R) , < t < T . 
We put A(a) := l( ao .oo)( a ) e<1 ^ r so that 

/>0O 

A u (t,x) = / A(a) u(t, a, x) da , (t,i)e[0,T]x(l, 



(23) 



Then we deduce from (JTSJ), (HU), (JUD, OH), and ((231) that, for < s < t < T < 1, 

l|A-*(«)(*) - A *(u)( s )llw 2 ' 5 



< 



I^A„(i,a)|| £ (w^ B ) ll£(«u(a)K(a)||w™ e( * a)/rda 



- So " H(r)dr _ - f g s " n(r)dr 



\\UA a (t,a)\\ c{W 2^ U(v u (a))v u (a)\\ W M e (t a)/r da 



||f/ A „(t,a) - t/A tl (s,a)|| £ (^ B ,H/ p M B ) ll£K(a)R(a)llw^ e (t a)/r da 
+ / WA u (s,a)\\^ W 2^\\^(v u (a))v u (a)\\ W 2^\X(t - a) - X(s - a)\da 

•J 

/» oo 

+ \\U Au (t,0) - U An (s,0)\\c(w*« B ,w™) I e (o+t '/ T da 

+ II^a„(s,o)||£ (w ** 

+ ||C/a u (s,0)|| £ ( W t2^) / ||u°(o)||w|tf |A(a + i) — A(a + s)| da 

« 

< c(R) (t-s) + c(R) (t - sy-o + c(R) [ \X(t - a) - A(s -a)\da 

Jo 

/•OO 

+ c(R) / \\u°(a)\\ W 2i> |A(a + t) - A(a + s)| da . 



|u°(a)llw™ e (a+t)/r da 



Next note that 

|u°(a)|| H A2a |A(a + t) — A(a + s)| da 

"(oo-s)+ 



< 



|«°(a)|| W 2* e^+^da + 



(a -t) + 



|«°(a)|| 



U" 2r ' 



(a -«) + 



e (a+t)/r _ e (a+s)/r 



da 



< e^)/r \\ u °\\ Lgmao) , wr) |(ao - »)+ - (ao - t)+| W + r~ 2 e 1 /- (i - s ) \\ u °\\ 

< c(R)((t- s y/"' + (t-s)) 
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owing to p6p while 

/■(s-oo)+ p(t-a )+ 

\\{t-a) -A(a-o)|do < / ( e (*- a )/^ - e (— «)/^) da + / e {t ~ a)/T da 

JO J(s-a ) + 

< r- 2 e 1 ^ (t - s) + e x ' T ((t - a ) + - (s - a )+) 

< (1 + r" 2 ) e 1 ^ (t - s) . 

Therefore, 

\\A Hu) (t) - A Hu) (s)\\ W 2, <\t-s\ K , Q<t,s<T, 

due to the choice of k provided that T = T(R) G (0, 1) is chosen sufficiently small. Furthermore, 
using (H]), (HU), and $Z3§ we obtain for < t < T 



\Mu)(t)\\ Ea < / \\U Au {t:t~a)\\ c(W 2 :B) \\i{v u {t-a))v u (t~a)\\ wr e a / T Aa 

J 

/"OC 

+ / \\UAAt,0)\\ c{w 2 S W 2„ } \\u a (a - t)\\ w 2 B e a / T da 



< c{R)T + (c(R)t^ +c )ey T \\u a \\ Ee 

< 1 + R 

provided that T = T(R) G (0, 1) is chosen sufficiently small. Since is obviously non-negative 
and <I>(ii) G C([0, T], E a ) holds by similar arguments as used to prove the Holder continuity of A$( M j, 
we conclude that $ maps Vt into itself. That it is a contraction follows from the observation that 
if u, u G Vt and < t < T, then 

||*(u)(t)-«l>(ii)(*)lk 

< / l|t/A„(*,*-o) - C^ a (*,*-a)||/:(^ B ) U(v u (t - a)) v u (t - a)\\ W 2„ e a/r da 

+ / WAu(tit- a)\\c(w 2 %) U( v u(t - a))v u (t - a) - £(u fi (t - a))v a {t - a)\\ W 2„ e a/r da 
Jo P ' B p 

/OO 
||^ u (*,0) - U A ,(t,0)\\c { w^,w^) W u °( a - *)llw^ B e<l/Tda 

and hence, using (TTJ]), (13]), ([11]), (JTHD , (H, and ((231), 

/"* 1 

||$(u)(t) -$(u)(t)|| E<r < c(R) \\u- u\\vr / e a/T da + c{R)t - a \\u-u\\ Vt < -\\u-u\\ Vt 

Jo * 

provided that T = T(-R) £ (0, 1) is chosen sufficiently small. Therefore, by Banach's fixed point 
theorem there exists a unique u G Vt such that $(u) = w. Note that (TT8)) and (|23p imply 
u = $(u) G L oo ([0,T],£:i) since u° G Si, whence u u € C([0, T], VF p 2 e ) by Lemma O Due to 
u(T) G v u (T) G W 2 g , and the fact that T was chosen depending only on R satisfying (fT6|) . we 
can iterate this argument and extend u and v u uniquely to functions u G C( J, E^) PI L^^^J, Ex) 
and u G C^J, Wfy) n C( J, I4^ 2 e ), where t+ := sup J = oo if 

sup {|K«)|k + ||«(t)|U e ((o,ao),w?')} < f or a11 T>0. (24) 

0<t<min {t+,T} 

Clearly, this so extended function u still satisfies 

/ e~^^^U A (t,t-a)£,{v(t-a))v{t-a) , a < i , 
~ 1 e-/-^ d ''^(t,0) U °(a-i) , a>i (25) 
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for a > and < t < t + , where we simply write A — A u and v — v u . Next recall that 
u° G E x n C a (R + ,L p ) and so, for t G (0, t+) and a > with a ^ t, 

d t u(t, a) = l [a<t] (t, a) e- # "W dr { - A(t) U A (t, t - a) £(v(t - a)) «(t - a) 

+ U A (t,t-a)(d t + A(t-a))(Z(v(t-a))v{t-a))} 
-l [o>t] (t,a) [fi(a-t)u(t,a) + e~ ^-^ (r)dr A(t)U A (t,0)u°(a - t) 
+ e - /.% ^ dr U A (t, 0) d a u°(a -t)\ 

and 

d a u(t, a) = — /i(o) it(t, a) 

- l [o<t] (t, a) e~ fo MMdr ^ t _ a) ( A(t _ fl) + ^ _ v{t _ a) ) 

+ l [o>t] (f,a) {e-^-^ Wdr C/A(t,0)a aU °( a -t) + u(a-*)n(*,a)} . 

Thus is a solution to (H])-© with the regularity properties as stated in the assertion of the 

theorem. 

It remains to prove that t + = oo. Wc fix T > arbitrarily and put Jt '■= Jn [0, T], Defining 
A4i(t, x) :— J °° e a l T u(t, a, x) da we observe that © and (jTO]) ensure 

d t v(t) < v(t) + \\fi\\ooMi{t) , teJ, 

whence 



J ||A^i(s)||cods + . i e Jt 



ll«(*)l|oo < c(T) 
But since 

A4i(t) < ||ei|ooe T/T / ||«(t -a) Hoc da + c(T) ||t»°|| Bl , t G Jt 



by (|22|) , we conclude 

||«(t)||oo + ||A4i(t)||oo < c(T) , i G Jt ■ (26) 
Owing to u° G C([0, a ], V^- 1 ) and u G C( J W£ B ) it follows from 

\\U A (t,s)-U A (t,r)\\ c(w2 w 2^-i ) <c(t )(s-r) 3 ^ , < r < s < t < t < t+ , (27) 

and ((25l) that u(-,a ) G C( J \ {a }, T4^ 2w_1 ). Property (J27J) is shown analogously to [4j Equa- 
tion (II. 5. 3. 8)]. Provided a < f+, (JTSJl warrants that lim t/Q0 u(t, a Q ) = u°(0) in W p 2l/ with 
n/p < 2z/ < 2cj - 1, while {27]) warrants that lim t ^ ao u(t, a ) = £(v°)v° in Wg"' 1 . Thus, the im- 
posed compatibility condition on u° and u° entails u(-,a ) G C(J, W^") c — > C(J, C(f2)). Recalling 
that A u is given by ([6]) we set 

1 f 00 

/(i,x) := e ao/T u(t, a ,x) H A u (t,x) - / e a/r u(a) u(t, a, sc) da , (t, x) G Jt x 17 , 

and deduce / G C(Jt X fi) with 

|/(i,x)| < c(T) , ((,i)eJ T x!i (28) 

due to (f2T>l) and 

N«.*)i-s{is?i-_'sr:" D)l "' ^:}^< r »- 

We then observe that A = A„ solves the quasilinear parabolic problem 

d t A - div x (£)(A)V x A) = f(t,x) , (i,a:) G J T x f2 
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subject to d v k = and A(0) G W^ B . We refer to ©, fj26|>. and ((2Sj) when using Lemma 5.1(ii)] 
to obtain that A G BUC s (J T ,C 5 (0,)) for some <5 > 0, and hence A G BUC £ (Jt, C 1 (f2)) for some 
e > by [51 Lemma 4.2, Remark 4.3], where BUC £ stands for 'bounded and uniformly e- Holder 
continuous'. But then A = A u is uniformly Holder continuous, that is, 

\\A u (t) - A u {s)\\ c{W 2 B>Lp) < c(T)\t~s\ £ , MG Jt , 
so that [H Lemma II. 5. 1.3] implies 

II^A u (M)IW eR) + (t-sy\\U Au (t,s)\\ 2e < c(T) (29) 

\ p ,B 1 \ P ' p .13 s 

for t,s G Jt with s < t. Note that c(T) depends here on T only (but not on some norm of u). 
Combining (gU), (ESI), and dm we have 

/■mill {i .ao} 

ll M WIH 9 ((o,ao),w-) < / W U A(t,t-a)\\ g ciLp>w2v U(v(t-a))v(t-a)\\l p da 



\\U A (t,0)\\l {w * %) ||« (o-*)[[^« da 

min{t,ao} p ' P 



/•min {i,ao} 

< c(T) / a-^da + c^H^lli^^^),^) 

Jo 

< c(T) 

for t G Jt thanks to < 1. Finally, from (|2"5)) . ([2"u]) . and it follows analogously that 
it G Loo(Jt, Eg). From this and (j2"4"|) we deduce i + = oo. This proves the theorem. □ 

3. The Degenerate Case 

We now turn to the "degenerate" case where D is allowed to vanish but only for A = 0. More 
precisely we assume that D G C 2 ~(M) is such that D(0) = 0, D(z) > if z > 0, and 

I D ~ f 1 Z D '}f d z < oo and lim Z ! f 'ffj = . (30) 

Jo D(z) z^o D(z) 1 / 2 1 ' 

The function 

$Z5(z):^2; / — - dy- / . dy , z G [0, oo) , (31) 



D(y) » J Q D(y) 

is then a well-defined smooth convex function satisfying $d{z) > <&£)(!) = — Id for z G [0, oo). We 
also put 

D(z) := / D(y)di/ and := / S(i/)dy, z G [0, oo) . 

Jo Jo 

We note that both D(z) = z m ~ x , m > 1 and £>(z) = e" 1 / 2 fulfil (|30j) . 

As for the differentiation and dedifferentiation rates £ and /j, we assume that £ fulfils ([5]) while 

fi satisfies 

fi G W^,(R) is non-negative and n(a) = for a < a , (32) 
the latter assumption being stronger than (fT0|) . 

Finally, the initial data are required to satisfy the following properties: 

< u° G Li((0,oo) x n,e a/T dadx) n £«,(((), oo) x fl) and < u° G L OQ (Q) (33) 

and 

A G £oo(tt) and D (A ) G W^Q) , (34) 

where 

poo 

A°(x) := / e a ' T u°{a,x)da , (35) 
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Theorem 3.1. Let T > and put U := (0, T) x (0, oo) x fi. There are two non-negative functions 
u£L 00 (U) andu e Loo((0,T) X Q) nC([0,T],L 2 (n)) satisfying 



1 f°° 
d t v(t,x) = -(1 - a;))) u(i,a;) + / e a/r /x(a) u(t, a, x)da a.e. in (0, T) x O , 

and 



T 



= / / / a, x) (<9tV> + 9 a -0) (i, a, x) dadxdt + / / u°(a, x) ip(0, a, x) dadx 

Jo Jn Jo Jn Jo 



1 <- T 



T 

r-T 



t;{v(t, x)) v(t, x) ^{t, 0, x) dxdt — Iff fj,(a) u(t, a, x) tf>(t, a, x) dadxdt 
o Jn Jo JnJo 



pi p p oo 

/ / / J(<, a, x) \7 x i/j(t, a, x) dadxdt 
Jo Jn Jo 



for if) £ C^([0, T) x [0, oo) x 0) mt/i d„V(*> a,x)=0 for (t, a, x) £ (0, T) x (0, oo) x dfl, where the 
functions A and J are given by 

poo 

A(i,x) := / e a/T u{t , a, x) da for a.e. (t, x) £ (0, T) x SI , (36) 



J := V x (uD(A)) - uV x D(A) in V'(U,R n ) , (37) 

and satisfy A £ C([0,T], L 2 (Q)), D(A) £ L 2 ((0, T), W^(Q)), and J £ L 2 (U,R n ). 

The proof of Theorem 13. II is performed by a compactness argument, approximating the diffu- 
sivity D by non-degenerate diffusivities (D a ) a> Q for which we can apply Theorem 1 2 . 21 and obtain 
a sequence of solutions (v a , u a ) a> Q. The next step is to pass to the limit as a — > and we now 
point out the difficulties to be overcome: first, as D(A) vanishes when A = 0, the equation |T]) 
is no longer uniformly parabolic with respect to the space variable and V x u is unlikely to be a 
function. Furthermore, as a > 0, we may have A(i, x) = but u{t, a,x) ^ for a £ (0,ao) and 
(|T|) gives no information on u in that case. We therefore cannot expect to have strong convergence 
on the sequence (u a ). There are however nonlinear terms in |T]) and ^ for which strong conver- 
gence is necessary to identify the limit. In particular, the strong compactness of (v a ) a is needed 
to pass to the limit in the term £(v Q ) v a . As v a solves an ordinary differential equation, such a 
compactness can only be obtained as a consequence of that of (t, x) i — > J °° e°/ r /i(a) u a (t, a, x) da. 
One step in the proof is thus to show that certain integrals of u a with respect to age enjoy some 
compactness properties with respect to the time and space variables. The strong compactness of 
A a (t, x) = J°° e°/ r u a (t, a, x) da will also follow from this step. Next, in order to identify J, strong 
compactness is needed on (\J x D(A a )) a to pass to the limit in the term u a V x D(A a ). This is proved 
by a suitable adaptation of an argument from [5] . 

We now begin the proof of Theorem 13. II We fix T > and consider a sequence (D a ) a £(o,i) of 
functions in C 2 ~(R) with the following properties: for every a £ (0, 1), there is d a > such that 
D a {z) > d a for all zGl, and 

D a (z) = D(z) for ze[ae" T|WI °=,co). (38) 

Next, let (w° , u° a ) Qg ( ,i) be a sequence of non-negative initial data fulfilling all the requirements of 
Theorem 12.21 together with the following properties: 



lim|y J e a/T \u° a (a,x) -u°(o,a)| dadx + \\v° a -v ]]^ =0 , (39) 

and there is cq > such that 

/ / e Q / T <(a,x)dadx+||<|| oo + ||t;°|| oo + ||A^|| oo + ||i3(A°)|| lvl <c (40) 
J si Jo ' ' " 2 ' 
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with 

t>oo 

A° a (x) := / e a ^u° a (a,x)da>a (41) 

Ja 

for all x G Vt and a G (0, 1). 

We denote by (v ai u a ) the solution to 

d t u a +d a u a = div x (£) Q (A Q ) V x m q ) - fi(a)u a , (t, a, x) £ (0, oo) x (0, oo) x , (42) 

1 f°° 
d t v a = - (l - £(w a )) v a + / e a/T /x(a) u a (t, a, x) da , (t, x) € (0, oo) X Q , (43) 



where 

A a (t,x) := I e a/T u a {t, a,x)da , (i, x) 6 (0, oo) X Q , (44) 



subject to the boundary conditions 

u a {t,0,x) = -£(v a (t,x))v a (t,x) , (t, x) £ (0, oo) x O , (45) 

T 

£»a(A Q ) 9* tt Q = , (t, x) e (0, oo) x ffi , (46) 

and the initial conditions 

it a (0, a,x) = u°(a,x) , w Q (0, x) = u° (x) , (a, x) G (0, oo) x O . (47) 
We note that, thanks to (|42|) . A a solves 

d t A Q = div x (Z> a (A a )V x A a ) +g 1 a - g 2 a in (0,T)xfl (48) 
with homogeneous Neumann boundary conditions and 

A ft xl f 00 
g^(^x) := e ao/r u a (t,a ,x) + ° V ' y >0 and g*(t,x) := / /i(o) e a/r u a (i, a, x) da > . 

For further use, we introduce the following functions: 

/>oo />OC 

M l a (t, x) := / e a / r u a (i, a, x) da and M.2 a (t,x):= / e a ^ T /i(a) u a (t, a, x) da 



for (i,x) G (0,T) x Q and a £ (0, 1). Using again (j42|) together with (J44J) , we realize that Mi, a 
and A^2,a are solutions to 

9tMi, a = div x (D a (A a )V x Mi, a ) +-£(v a )v a + — — - / n(a) e a/T u a {., a, .) da (49) 



and 



d t M 2 , a = div x (D a (A a )V x M 2 , a ) + + / (//(a) - ^(a) 2 ) e a/r u a (., a, .) da (50) 



in (0,T) x J7, respectively, with homogeneous Neumann boundary conditions. 

In the following, we denote by c and Cj, i > 1, positive constants depending on D, fi, £, ao, 
r, and Co in (|40p . but not on a. The dependence upon additional variables (such as T) will be 
indicated explicitly. As in the non-degenerate case, we establish Loo-bounds for A4i, Q , v a , and u a . 

Lemma 3.2. For a £ (0, 1) and t £ [0, T], we have 

Wlloo < cr(T) , (51) 

pT p poo 

/ / / D a (A a )(t,x)\V x u a (t,a,x)\ 2 dadxdt < c x (T) . (52) 
Jo Jo JO 
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Proof. On the one hand, since fj, belongs to Loo(0, 00), we have 

& > llMlloc 

T 

by (pL3"j) . from which we deduce that, for t 6 [0,T], 

v a (t,x) <vl{x) e^ + IHloo / Xi. Q ( S ,x) e ( t - s '/ T d S . 

Jo 

Using ([40|) gives 

||«o(*)||oo < c(T) (^ 1 + J H-^l 

,a (^) || oo 

ds^j . (53) 
On the other hand, it follows from (pf5|) and the boundcdness of £ that 

d t M lta < di V:c (Z> a (A a )V x A4i >a ) + + ^ . 

The comparison principle then ensures that 

Mi, a (t,x) < \\M lia {0)\\°o e t/T + l J ll^(s)||ooe(*- s )/ T d S < c(T) ^1 + ^ KOOIUd^ 

for t G [0, T]. We now combine this estimate with (f53|) and end up with 

||Ml,«(i)||oo < c(T) ( 1 + J Q J Q ll-Mi^^Hoodad^ < c(T) ^1 + ||X ljQ ( s ) ^ d^ 

for t E [0, T]. The Gronwall lemma then gives the claimed bound on ||Ati iQ ||oo, which in turn gives 
that for \\vJoo by 

Finally, by l[42j). ([45]) . and l|i6]). u Q satisfies 

9 t w« + <9 a u Q < div x (-D a (A a ) V x u a ) 

with it a (i, 0,x) = t;(y a (t,x))v a (t,x) j r and subject to homogeneous Neumann boundary condi- 
tions for t € [0, T]. On the one hand, the comparison principle readily implies that 

u a (t,a,x) < \\u° \\oo + - sup Hua^lloo , (t, a, x) e (0, T) x (0, oo) x fi , 
T se[o,T] 

which, together with and the already established bound on ||u a ||txj, allows us to complete the 
proof of (|5ip . On the other hand, since u a is non-negative, we also have 

/> />oo pT /• poo 

/ / u a {T) 2 dadx + 2 / / / D a (A a (t)) \V x u a {t)\ 2 dadxdt 
Jn Jo Jo Jn Jo 

< (u° a ) 2 dadx+ / / u a (tf dadxdt + / (£(«„(*)) v a (t)) 2 dxdt 
Jn Jo Jo Jn Jo Jo Jn 

< iKlloo Klli + / T |Mt)||oolK(«)||idt+ [ T \\v a (t)\\lM 

Jo Jo 

< c(T), 

the last inequality being a consequence of (|51[) . and the obvious inequality ||u a (i)||i < 

||-Mi,a(*)||i forte [0,3*1. □ 

We next derive some estimates for A Q . 
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Lemma 3.3. For a G (0, 1), t G [0,T], and x 6 Q, we have 

l|Aa(*)IU+ / [ \v x D(A a )\ 2 dxds+ [ [ \ \ ■ - ,-._,( :n . (.-,4) 

Jo Jn Jo Jn 



and A a (t,x) > ae" T|W ". 



n Aa 

D(A a )(t)\\ lro .+ [ T \\d t D(A a )( s )\\lds < c 2 (T), (55) 







A straightforward consequence of (|38|) and the last assertion of Lemma 13.31 is that 
D a (A a )(t,x) = D(A a ){t,x) for (t,x) e (0,T) x fi . 

Proof. Clearly A Q < A4i, Q and the Loo-bound for A Q is a straightforward consequence of Lemma l3.2 
It next follows from [|4"5]) that <9iA Q > div x (L> a (A Q ) V X A Q ) — ||/i[[oo A a in (0, T) x Q with homoge- 
neous Neumann boundary conditions. As t ^ a e~* is a subsolution to the previous equation, 
the lower bound A Q > ae _T " M "°° in (0, T) x SI readily follows from pTj) by the comparison prin- 
ciple. 

We next multiply |@5} by <J>' D (A Q ) with being defined in (|3"Tj) and integrate over (0,T) x f2 
to obtain 

f ($ D (A a (t))-$ D (A° a )) dx+ [ [ <S> D (A a )D(A Q ) \V x A a \ 2 dxdt 
Jn Jo Jn 

= f [ <7a $ij(A«) dxdt - / f g 2 a $' D (A a ) dxdt . 
Jo Jn Jo Jn 



On the one hand, since $>' D is non-positive in (0, 1) and g\ > 0, we infer from (|5ip and the 
Loo-estimate on A a that 

f [ ffi^CAaJdardf < [ [ gll [ll00) (A a )& D (A a )dxdt 
Jo Jn Jo Jn 

' e aalT \\u a \\oo + Jl^W + ||A Q ||oo) dxdt 



<- S\ 

Jo Jn 
< c(T) 



On the other hand, as & D > on (1, oo), & D < on (0, 1), and r \& D {r)\ < Id for r G [0, 1], we 
have 

pT p pT p poo 

gl& D {A a )dxdt < n{a)e a/T u a (t,a,x) dal (0 ,i)(A Q ) \& D (A a )\ dxdt 

Jo Jn Jo JnJao 

< \\fi\\oo [ [ Aa,l[o,i)(Aa) |^£>(A )| dxdt 

Jo Jn 

< WnW^ I D T\n\ . 

Recalling that $u(r) > — Id for r > and D& D = (D') 2 , we conclude that 

/ f \\7 x D(A a )\ 2 dxdt < c(T) + f $ D (A a (0)) da; < c(T) + max {I D , $n (c )} • 
Jo Jn Jn 
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Similarly, we multiply (|4"5)) by log A a and integrate over (0, T) x ft: using the non-negativity of 
g\ and , (|4T))) . and the Loo-bound on A a we obtain 



/ / ^ |V.A a | 2 dxdt 
Jo Jn A a 

< [ A° a (logA° - 1) dx - / A a (T) (logA a (T) - l) dx 
Jn Jn 



in 



+ [ [ g x a logA Q l [ i )0o) (A Q )da;dt- / / g 2 a logA a l( ,i)(A Q ) dxdt 
Jo Jn Jo Jn 

< c(T) + T |0| Halloo log (1 + || A Q ||oo) + ||a*||oo / / A a |logA a |l (0 ,i)(A a )da;dt 

Jo Jn 

< c(T) 

as ||5q||oo is bounded uniformly with respect to a £ (0, 1) by (|51[) and the Loo-bound on A a . 

We next multiply (gHJ) by 2d t D(A a ) and integrate over (0,t) x fl, t £ [0,T]: using (gDJ), (jHj) , 
and (l54f we obtain 

2 / f D(A a )\d t A a \ 2 dxds + \\V x D(A a ){t)\\ 2 2 
Jo Jn 



li ||V,DiA:;)|!;; +2 / / D(A a ) |9t A Q ^e o °/ T u Q (3,a ,a;)+ f - + M« | A„ j dr,b 



o Jn 



< c(T)(l 



< 

Therefore 



/ / D(A Q ) |9jA | dxdi 
Jo Jo 

/ / D(A a ) \dtA a \ 2 dxds + c(T) . 
Jo Jn 

f f D{A a ) \d t A a | 2 dxds + || V, D(A a )(t) || I <c(T), 
Jo Jn 



from which the claim (55) follows as |d t D(A Q )| < c(T) y/D(A a ) \d t A a \ by (53). □ 

At this point we have gathered the information required to show the strong compactness of 
(A a ) Q . This is, however, not sufficient to pass to the limit as a — > as there is a nonlinear 
dependence on v a in (|43[) . We now aim at proving the strong compactness of (v a ): this will be 
achieved by the strong compactness of (M 2 ,a)a which we show now. 

Lemma 3.4. For a £ (0,1), t £ [0, T], and S £ (0, 1), we have 

(\\Vx(M 2 , a -6)l\\l+\\d t (M 2 , a -5) 2 + \\ wi+i(ny ') dt<c 3 (T,S) . (56) 
Proof. We multiply ([50]) by (.M2, a — S), and integrate over (0, T) x 17 to obtain 



o 



1 ' T 



r || (^2,a(T) - 5). L + / / D(A a ) (V, (M 2 , Q - <S). | dzdi 
^ Jo Jn 

< ±\\(M 2 , a (0)-6) + \\l + T\n\ ( l| - M2 T a||oo + ||/x / || 00 ||Mi,a|| c) \\(M 2 . a -5) + \\ oo 
As (A^2,a — #) . < M-2,u < IImIIoo A4i !a and fi £ W^,(0, oo), we infer from Lemma I3T21 that 



/ / D(A a )|V B (A4 a ,a-*) + | 2 daxi* <c(T) . (57) 
Jo Jn 
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Now, on the one hand, since the support of [i is included in [do, oo), we have M. 2 ^ a < \\[i\\oc A a 
and 

{(t,x) G (0,T) x fi : M 2 , a (t,x) >5}d {(t,x) e (0,T) x : A a (t,x) > <5/IHIoc} • 

Introducing mg := min^/n^n^ oo) D > we deduce from (|54[) , (|57[) . and the previous observation 
that 



/ ||Vx(A<2,o-<5)+ ||ad* = 4/ / (A^ 2 , Q ™<5)+ | V, (X 2 , Q - <5) + | 2 dxdt 
Jo Jo Jn 





— r 


< 






m s Jo J 


< 


c(T,5) [ 




Jo 


< 




in 





Hoc Aa - *)+ D{A a ) I V* (M 2 , a - S), I dxdt 



I I D(A a )\V x (M 2 , a -S) + \ 2 dxdt 
Jo Ja 



On the other hand, if ip € W^ +1 (fl), it follows from ((50)) and Lemma [3721 that 
9t (A^2,q - 5) 2 + tpdx 

{M 2 , a - S) + d t (M 2 , a - S) + Tpdx 

< 2 f {M 2 , a - 6) |V S VI ^(A a ) |V :c X 2 , Q | dx 
Jn 

+ 2 / \iP\D(A a )\V x (M 2 , a ~S) + \ 2 dx 
i>h\\(M 2 , a -5) + tt 



-||M 2 ,a| 
T 



+ 2 (Halloo + IHlL) ll^i.alU Will || (M 2 , a - 5) + IL 

< 2||A<2,a|L l|V^||2p(A Q )||^ 2 ^D(A a )\V x (M 2 , a -5) + \ 2 dx 
+ 2|H| M / Z>(A Q ) IV,, (7W2, a - 5) + | 2 d2T + c(T) H^lli 

< c(T) (||V^|| 2 + Halloo) + jf -D(A Q ) |V x (A4 2 ,a - S) + \ 2 dx 
Owing to the continuous embedding of W* +1 (f2) in L 00(H) we conclude that 

\\d t (M 2 , a - 5)1 \\ w t +i{ny < c(T) (l + D(A a ) \V X (M 2 , a - S) + \ 2 dxj 
which together with (|57|) implies the second claim in |5 



1/2 



□ 



Lemma 13.41 provides the desired compactness for (M.2,a)a with the help of the following lemma. 

Lemma 3.5. Let Q be an open bounded subset ofVi N for some N > 1 and p 6 [1, 00). We consider 
a sequence (zk)k>i of non-negative functions in L p (Q) and assume that there is a sequence (Zj)j>\ 
in L p (Q) such that 



lim 







HI 









= for all j > 1 



(58) 
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Then (zk) converges in L p (Q) as fc — > oo. 
Proof. For i > 1, j ' > 1, and fc > 1 we have 



I^-^IIp < 



< 



Zk - - 

I 



Zk - - 
I 



Zk 



101 



1/p 



Zk r 

3 



Zk r 



Z, 



Letting fc — > oo and using 1|58|) give 



\Zi-ZjWp < IQI 



i/p 



so that (Zj) is a Cauchy sequence in L P (Q) and there is Z <E L p (Q) such that 



lim — Z 



Next, for j > 1 and fc > 1, we have 



\zk - Z\ 



< 



< 



101 
i 



. 



Zk - 



(59) 



\Zj-Z\ 



Zk - - 
J 



-z, 



\Zj-Z\\ p 



hence 



101 



\z 4 



limsup \\z k - Z\\ p < , _ , ; , 

/c — >oo J 

by (fS"5)) . Letting j — » oo and using (fS"9"|) give the expected convergence. 



□ 



Finally, to link the limits of (A a ) a , {M.\ t a)a and {M.2,a)a with that of (u a ) a we need to control 
the behavior of u a for large a and report the following result in that direction. 



Lemma 3.6. For a G (0, 1), t G [0,T], and 4> 1 idc /lave 

e a / r u Q (f, a, x) dadx < Cs(T) u) a (A) with uj a (A) := 



Q J A 



1 



n ja/2 



e a/r u„(o,x) dadx + — 



Proof. Let 77 G C°°(R) be a fixed non-decreasing function such that 77(a) = for a < 1/2 and 
77(a) = 1 for a > 1. For A > 1, we multiply (j4l?l) by r](a/A) e°/ T and integrate over (0, 00) x Q with 
the help of (|4"6")l . Since 77(0) = we thus obtain 



n Jo 



'/ 



e a ^ T u Q (i, a, x) dadx , 



-t/r 



n Jo 



77^—^ e a / r u a (t, a, x) dadx^j < ~~~^ 2E J Mi,a{t,x)dx 



By virtue of (|51[) the right-hand side of the above differential inequality is bounded by c(T) /A and 
Lemma 13.61 follows after time integration, taking into account the properties of 77. □ 

Proof of Theorem^ Recall that (D(A a )) a is bounded in L oo ((0, T), (ft)) and (<9 t Z)(A Q )) Q 
is bounded in L 2 ((0,T) x fi) by |55|) . Owing to the compactness of the embedding of W% (Q) in 
L2 (f2) we may apply [15l Corollary 4] to conclude that 

(D(A a )) a is relatively compact in C([0,T],L 2 (Q)) . (60) 
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A similar argument allows us to deduce from Lemma EUand [13 Corollary 4] that ((M 2 , a - Vj)+) c 
is relatively compact in L 2 ((0, T) x Q) for each j > 1. Since (A^2,ce)„ is bounded in L oo ((0, T) x Q) 
by (f5"Tj) , the Lebesgue dominated convergence theorem actually allows us to conclude that 

((M2,a — l/j) + ) a is relatively compact in L2((0,T)xfi) (61) 

for each j > 1. We then infer from Lemma [3T2l ((54)) , and (f60|) that there are a sequence (ak)k>i, 
a k -> 0, three functions £ G C([0, T], £ 2 (0)) n L 2 ((0, T), W 2 \n)), d G L 2 ((0, T), W^fi)), and 
u 6 Loo(C/), and a sequence (Wj-)i>i m -^((O, T) x f2) such that 



(w Q J fe A u in ioo([/) , 

(£>(A a J,S(A a J) fc -(d^) in L 2 ((0,T),Wl(n)) 

(D(A ak )) k ^£ in C([0,T],L 2 (fi)) , 

((A4 2>afc - l/j) + ) k — » W, in L 2 ((0,T)xO). 



(62) 
(63) 
(64) 



Combining the last convergence and Lemma [3.51 actually give that there is W <E L 2 ((0, T) x f2) 
such that 



{M %ah ) k ^W in L a ((0,T)xfi) 



(65) 



In addition, as the function I? is a diffeomorphism from (0, oo) onto its range with inverse D , 
the bound ([M]) and the convergence imply that 



(A a J k — IT V) m C([0,T],L 2 (f])) 



(66) 



We now claim that, if x is a non-negative measurable function such that x( a ) Se a / T for a.e. 
a > and some 3 > 0, we have 



with 



M x , a (t,x) 



(M x , ak ) k ^M x in M(0,7)xO) 



x(a) u a (t, a, x) da and .M x (i,x) 



(67) 



x(a) u(t, a, x) da 



o 



for (t,x) G (0,T) x fi. Indeed, consider t/> G Loo((0,T) x fl). For A > we have 

{M Xt a k - M x ) (t, x) ij}{t, x) dxdt 

I II (u ak — u) (t, a, x) rp(t, x) x(a) dadxdt 

(u ak — u) (t, a, x) tp(t, x) x(a) dadxdt 



o Jn Jo 

T 



< 



o Jn J a 

If 

o Jo. Jo 



(u ak — u) (t, a, a;) ip(t, x) x(a) dadxdt 



s a / T (u ak + u) (t, a, x) dadxdt 



lo JaJA 

We then infer from (|39| , Lemma 13.61 and (|62|) by a weak convergence argument that 



z a/r u(t, a, x) dadxdt < Tc 5 (T)u (A) with w (A) 



o Jn J a 



n J a/2 



e a/T u°(a,x) dadx+— . 
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It also follows from (|39[) and Lemma that 



[ [ e a/T u ak {t,a,x)dadxdt < T c 5 (T) ( I I e a/r \u° ak - u°\ dadx + uj (A) 
o J o J A \Jq J a 

and thus 

limsup / / / e a/T u ak (t,a,x)dadxdt <Tc 5 (T)lo q (A) . 

k^oo Jo Jn J A 



Consequently, by (|62|) . 



lim sup 

k — >oc 



(M x ak — M x ) (t, x) tp(t, x) dxdt 



o Jn 



< C (T)S||V||ooWo(^) . 



Since the above inequality holds true for all A > and tuo(A) — > as A — > oo by (|33|) . we may let 
A — > oo to conclude that (-A4 XiQfc ) fc converges weakly towards Ai x in Li((0,T) x f2) as fc — > oo. 
As < M Xtak < EMi, ah , (Mi, ah ) k is bounded in Loo((0,T) x fi) by l[5T]). and since (0,T) x £1 
has finite measure, the previous Li-weak convergence implies the claim (p7|) . 

In particular, we deduce from (|65p (with x( a ) = l[a .oo) ( a ) eQ ^ r an d x( a ) — eQ ^ T M( a )j re ~ 
spectively) that (A ak ) fc and {M.2,a k ) k converge weakly-* towards A and Ai 2 in £oo((0, T) x fi), 
respectively, with A and A4 2 given by 

/>oo />OC 

A(t,x):= / e a / r u(i, a, x) da and A4 2 (t, a;) := / e a ^ T fi(a) u(t, a, x) da 



for (t, x) £ (0, T) x fl. Combining this fact with ([65]) and (f66|) leads us to the identities A = D 1 (£) 
and J\4 2 = W, and we have thus shown that 

(AaJ fc — »A in C([0,T],L 2 (fi)) and (A4 2 , Q J fc — > A( 2 in L 2 ((0,T)xO). (68) 

A simple consequence of (|63|) and (|68|) is that d = £>(A) so that 

D(A)€L 2 ((0,T),W}(n)) and (£>(A Q J) fc - D(A) in L 2 ((0, T), W^fi)) . (69) 

We next denote by v the unique solution to 

d t v(t,x) = - (1 -£(v(t,x))) v(t,x) + M 2 (t,x) , (t,i)e(0,T)x$l, (70) 

T 

with initial condition v(Q) = v°. At this stage it is rather easy to deduce from (|43)) , (|70p. and the 
properties of £ that 

||K - < c 6 (T) (||« Qfc - «||| + ||M 2 , Qt - A4 2 ||i) . 

The strong convergence (JBSJ in L 2 ((0, T) x fi) of (A / ( 2 . Qfc ) fc towards A / ( 2 , ([M]), (gDJ), and the above 
differential inequality imply that 

(««Jk— * ^ C([0,r],Z 2 (fi)) . (71) 

Introducing 

J Q := D(A a )V x u a and j Q := D(A a ) 1/2 X7 x u a , 
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we infer from |@2J|, (gBJ, and |g7J that 

it a (i, a, x) (ftV' + c'oV') (i, a, a;) dadxdt + / / w° (a, x) ip(0, a, x) dadx 



o Jn Jo Jn Jo 



+ — I I £(v a (t, x)) v a (t, x) ?/>(i, 0, x) dxdt — III n(a)u a (t,a,x)ip(t,a,x)dadxdt 
T Jo Jn Jo Jn Jo 

pT /> poo 

J a (t,a,x)'Vxip(t,a,x)dadxdt (72) 

Jo Jn Jo 

for V G Cl([Q,T) x [0,oo) x Cl) satisfying d u ip{t,a,x) = for (t,a.x) G (0,T) x (0, oo) x 39.. By 
dSU) (ja)a is bounded in L 2 (?7, M n ) and so is J Q = £ ) (A Q ) 1 / 2 j Q by (|54"|) . We may then assume 
(after possibly extracting a further subsequence) that there is J 6 L2(U, K") such that 

(J Q J fe - J in i 2 (C/,M") . (73) 

Owing to (|39|) . (|62|) . ([7Tjl . and (|73| . we may pass to the limit as k — > oo in the weak formulation 
(1721) of gl]) to conclude that 

pT p poo p poo 

= / / / u(t, a, x) {ptij} + 3 a ip) (t, a, x) dadxdt + / / u°(a, x) -0(0, a, x) dadx 



o JnJo JnJo 



£(u(£, x)) v(t, x) tp(t, 0, x) dxdt — jjj fj.(a) ip(t, a, x) dadxdt 
o Jn Jo Jn Jo 



PI p POO 

/ / / J(t, a, x) S7 x ip(t, a, x) dadxdt 
Jo Jn Jo 



for V G Cl([0,T) x [0,oo) x 9) satisfying d v ip(t,a,x) = for {t,a,x) G (0, T) x (0, oo) x 50 as 
claimed in Theorem I3. II 

It remains to identify J: for that purpose we introduce the sets 

V := {(t,x) G (0,T) x 9 : A(t,x) > 0} , Z := {(t,x) G (0, T) x 57 : A(t,x) = 0} , 

and observe that J Q may be written 

J a = V x (u a D(A a ))-u a \7 x D(A a ) in P'((7,R") . (74) 

It follows at once from (|54[) . (JHSJ) , and the continuity of D that 

(u ak D(A a J) k ^uD(A) in L 2 (Z7) . (75) 

Next, we claim that, after possibly extracting a further subsequence (not relabeled), we have 

(y x D(A ak )) k — ► V X .D(A) in L 2 (C/,K") and a.e. in (0,T)xfi, (76) 

and adapt the proof of [5j Eq. (3.22)] to this end. We multiply flUJ) by -D(A a ) --D(A) and integrate 
over (0, T) x 9 to obtain 

f [ (D(A a )- 5(A)) 8 t A a dxdt = -/ [ \7 x (D(A a )~D(A)) ■ V x D(A a )dxdt 
Jo Jn Jo Jn 

+ f T [ (al ~ al) (D(A a ) - 5(A)) dxdt , 
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hence 

i-T 



f f V x (3(A a ) -5(A)) 2 dxdt = - f [ W x (D(A a )-D(A)) ■ W x D(A)dxdt 
Jo Jn Jo Jn 

+ f I (gi - gl) (D(A a ) - 3(A)) dxdt 

Jo Jn 

+ f f (5(A) - D(A a )) d t A a dxdt . (77) 
Jo Jn 

As V X 5(A) G L 2 ((0, T) x fi) by (J63J) and 5(A) = £ by ([66]) and ([68) we infer from (63]) that 

lim / / V x (D(A ak ) - 3(A)) ■ V X 3(A) dxdt = . 
fc ^°°Jo Jn 

It next follows from (51]) and (JMl) that ||^||oo + ||fi£lU < c(T). Therefore, by virtue of (Ml), ©, 
and (|68p we have 

lim / T / (ffi fc - S »J(B(A a J-D(A))dcd* = 0. 
k ^°°Jo Jn 

We next argue as in [5] (with D instead of z i— > z" 1 and X = W / 2 1 (^)) to show that 

lim / / (5(A) - .D(A a j) d t A ak dxdt = . 

Taking a = a& in (I77[) we may therefore pass to the limit as k — > oo and conclude that (T6"]) holds 
true. 

Now, as V x .D(A a J = (D'(A Qfc ) V x 3(A a k ))/D(A ak ) we deduce from (J68J) and ([76]) that 

^(A) 
D(A) 

Moreover, by ([54]) . 

|C'(A Q J|Ai{ 2 D(A ak )^ 

A l/2 

£'(A Q J| 2 A_ \ 1/2 ( f T f D(A ak ) ,„ A |2 V " 



(V x D(AaJ) fc — ► — V X 5(A) = V^(A) a.e. in V . (78) 



/ |V XJ D(A Qfc )| dxdi = / — — — — |V x A a J dardi 



2 A \ 1/2 



and the right-hand side of the above inequality converges to zero as k — > oo by ([30]) . ([68]) . the 
definition of Z, and the Lebesgue dominated convergence theorem. As V X D(A) — a.e. in Z 
by Stampacchia's theorem, we have shown that, after possibly extracting a further subsequence 
(not relabeled), (X7 x D(A ak )) k converges to \7 X D(A) a.e. in Z. Recalling ([78]) we have actually 
established that (V x D(A ak )) k converges to V X D(A) a.e. in (0, T) x SI which, together with (691) 
and the Vitali theorem, implies that (V x D(A ak )) k converges to V X D(A) in £i((0, T) x Sl,R ra ). 
Combining this convergence with (6"!?]) entails that 

(ua k S7 x D(A ak )) k — * wVifl(A) in Li(tf,R n ). (79) 

Thanks to (74]), (75]), and (79]), we conclude that J = W x (uD(A)) -uV x D(A) in V'(U,R n ) as 
claimed in ([37]) . In fact, as J and uV x D(A) both belong to L,2(U,M. n ), we realize that uD(A) 
belongs to L 2 ((0,T) x (0, oo), W%(Sl)). □ 
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Remark 3.7. It is actually not necessary to assume that D £ C 2 (R) and Theorem \3.1\ is valid 
if D £ C(R), so that it applies to the diffusivity 

D(A) = D max{A - A min , 0}" l_1 + e - 1/(eA) 

for m > 1 and e > 0. TTie proof is nevertheless slightly more technical as the sequence (D a ) a 
approximating D cannot coincide with D on some interval and has to be constructed carefully so 
that the above proof still works. 

4. Concluding Remarks 

A model accounting for the swarming of the bacteria Proteus mirabilis and involving age and 
spatial variables has been studied. It describes the evolution of small non-moving cells (swim- 
mers) and larger moving cells (swarmers), the latter moving according to Brownian movement 
divr (-D(A) V x uj with a diffusivity D(A) depending on the total motile swarmer cell biomass A 
defined in ([6]) and thus, in a nonlocal way (with respect to age), on u. Assuming that the dif- 
fusivity is bounded from below by a positive real number, existence and uniqueness of a strong 
solution has been established in Section [5] It is, however, expected on biological grounds that a 
certain amount of biomass is required for the motion of swarmers to be initiated, that is, D{A) is 
expected to vanish when A is below a threshold value A m ; n > 0. A step in that direction is made 
in Section [3] where the existence of a weak solution is obtained for A m ;„ = 0. To our knowledge, 
the more realistic case A m i n > has not been investigated so far, and we hope to return to this 
problem and to the formation of regular patterns as well in the near future. 

As a final comment, let us point out that in the model studied in this paper only Brownian 
motion is responsible for the movement of swarmer cells and describes somehow local displace- 
ments. Though, as pointed out in [9] and [5J, only swarmer cells of a certain maturity can actively 
participate in group migration, the so-called "raft building" , but nothing prevents young swarmers 
from being caught up in the flow and thus move with larger swarmers in the rafts. The diffusion 
term div x (D(A)\7 x u), however, reflects active movement of swarmers of any age, i.e. also of young 
swarmers. It is therefore more realistic to model migration by a drift term along the gradient of 
biomass, namely, d\v x (u E(A, v) V^A) with E > 0. The velocity E(A,v)V x A then points in the 
direction of increasing biomass density. The swarmer cell density equation including the above two 
spatial mechanisms then reads 

d t u + d a u = Aiv x (D(A) \7 x u + u E(A, v) \7 X A) - /j,(a) u , (t, a, x) £ (0, oo) x (0, oo) X fi , (80) 

instead of ([1]). The special case E(A,v) — D'(A) is actually stated in [6j. Note that the choice 
D = is possible in (|80|) and it would be interesting to see whether regular structures also arise 
from the model accounting only for drift motion. From a more theoretical viewpoint, the study of 
(ISP) seems to be more complicated than that of ([T]) because the initial-boundary value problem 
is no longer diagonal. Nevertheless existence of weak solutions can still be established and will 
appear elsewhere 
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